Impurity-Induced Quasiparticle Transport and Universal Limit Wiedemann-Franz 

Violation in d-Wave Superconductors 
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Due to the node structure of the gap in a d-wave superconductor, the presence of impurities 
generates a finite density of quasiparticle excitations at zero temperature. Since these impurity- 
induced quasiparticles are both generated and scattered by impurities, prior calculations indicate a 
universal limit (n 0,T 0) where the transport coefficients obtain scattering-independent values, 
depending only on the velocity anisotropy Vf/v2- We improve upon prior results, including the 
contributions of vertex corrections and Fermi liquid corrections in our calculations of universal limit 
electrical, thermal, and spin conductivity. We find that while vertex corrections modify electrical 
conductivity and Fermi liquid corrections renormalize both electrical and spin conductivity, only 
thermal conductivity maintains its universal value, independent of impurity scattering or Fermi 
liquid interactions. Hence, low temperature thermal conductivity measurements provide the most 
direct means of obtaining the velocity anisotropy for high Tc cuprate superconductors. 



PACS numbers; 74.25.Fy, 74.72.-h, 71.55.-i 
I. INTRODUCTION 

The characteristic feature of a d-wave superconduc- 
tor is the existence of four nodal points where the or- 
der parameter vanishes. Since low energy excitations are 
concentrated about these nodes, low temperature behav- 
ior is dominated by the details of the node structure 
and, in particular, the ratio of the Fermi velocity to 
the gap velocity (slope), Vfjvi- Prior theoretical work 
has shown that this velocity ratio is prominent in £x=. 
pressions for low temperature transport coefhcientsErEl 
as wellpas-the temperature dependence of the superfluid 
densitya"E3. However, discrepancies between values of 
Vflv2 obtained ham measurements of mic.ijnwave electri- 
cal conductivityO, thermal conductivityll3'E^, and super- 
fluid densityllil, as well as direct measurements of gap 
structure via ARPESEj indicate that the existing theo- 
retical predictions must be corrected through a more de- 
tailed analysis. To this end, we calculate herein electrical, 
thermal, and spin conductivity including the contribu- 
tions of vertex corrections and Fermi liquid corrections. 
Associated calculations of the superfluid density will be 
pursued in a future istuestigation. 

It has been shownE£l that for a superconductor with 
dx2_y2 pairing symmetry, the presence of impurities gen- 
erates a finite density of quasiparticle states down to zero 
energy (although the ultraJow energy regime remains the 
subject of some debatdlZnlS) . This results in a unique sit- 
uation where an increase in impurity density increases the 
density of quasiparticles while reducing the quasiparticle 
lifetime. As a result of the cancellation of these opposing 
effects, "bare bubble" conductivity calculations (neglect- 
ing the corrections we shall consider) indicate a univer- 
sal limit (O 0,T 0) where the transport coefficients 
attain constant values, independent of scattering.EI How- 
ever, we shall see that these results are modified by two 
types of corrections: vertex corrections and Fermi liquid 
corrections. Vertex corrections account for the fact that 



forward scattering does not interfere with the progress of 
a carrier to the same extent as back scattering. Hence, 
if the scattering potential varies in k-space such that the 
potential for forward scattering differs from that for back 
scattering, the bare bubble transport coefficients may be 
modified. Fermi liquid corrections account for the un- 
derlying Fermi liquid interactions between electrons in 
the superconductor. Due to such interactions, the pres- 
ence of a quasiparticle current induces an additional drag 
current which may renormalize the transport coefficients. 
The purpose of what follows is to improve upon the bare 
bubble results by including the effects of both types of 
corrections. 

In Sec. H, we define the parameters of our phenomeno- 
logical d-wave model, introduce the Green's function, and 
calculate the density of states. In Appendix neglect- 
ing all corrections, we derive a generalized bare bubble 
polarization function which can be applied to the cal- 
culation of either electrical, thermal, or spin conductiv- 
ity. By treating the general case, we avoid repeating the 
same basic calculation three times. In Appendix we 
calculate another generalized polarization function, now 
including the contributions of vertex corrections. The 
significance of vertex corrections in the universal limit is 
determined via a numerical calculation presented in Ap- 
pendix 1^. In Appendix we derive the renormalization 
of a generalized current due to the effects of underly- 
ing Fermi liquid interactions. In Sees. |l|, and 
we make use of the results in the appendices to calculate 
electrical, thermal, and spin conductivity in the universal 
limit (57 ^ 0,r ^ 0). Each of these sections begins with 
a derivation of the appropriate current density operator. 
These calculations reveal an extra gap velocity term in 
the thermal and spin currents due to the momentum de- 
pendence of the d-wave gap and therefore indicate « 
rection to the standard thermal conductivity formulalljH 
derived assuming an s-wave gap. Given each current op- 
erator, we present the bare bubble result and then note 
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modifications due to vertex corrections and Fermi liq- 
uid corrections. We find that contrary to the scattering- 
independent result obtained from the bare bubble calcu- 
lation, the universal limit electrical conductivity attains 
a vertex correction, which depends explicitly on the na- 
ture of the impurity scattering, as well as a Fermi liquid 
renormalization, which depends on the strength of Fermi 
liquid interactions. In addition, while the spin conductiv- 
ity is unaffected by vertex corrections (for small impurity 
density), it is renormalized due to Fermi liquid interac- 
tions. Only the thermal conductivity has neither a ver- 
tex correction nor a Fermi liquid correction. It therefore 
retains its simple, universal value. Conclusions are dis- 
cussed in Sec. ^ where we provide physical descriptions 
of the mathematical corrections calculated herein. 



II. D-WAVE MODEL, GREEN'S FUNCTION, 
AND DENSITY OF STATES 



To study the low temperature transport properties of 
a d- wave -superconductor, we employ a phenomenologi- 
cal modelld'l with the Brillouin zone of a two-dimensional 
square lattice (of lattice constant a), an electron disper- 
sion (via tight binding parametrization) 



efc = — 2ty (cos k^a + cos kyu) — ii, 



and an order parameter of dj, 
Ao, 



symmetry 



Afe = (cos k^a — cos kya) 



(2.1) 



(2.2) 



which crosses through zero at each of four nodal points on 
the Fermi surface, {k^ = ^ky). The key feature of such 
a model is that in the vicinity of each of the gap nodes, 
ek varies linearly across the Fermi surface and varies 
linearly along the Fermi surface. Defining local momen- 
tum variables at each of the nodes with ki perpendicular 
to the Fermi surface and k2 parallel to the Fermi surface, 
we can designate at each node both a Fermi velocity 



and a gap velocity 



V2 



Wyki 



Vf 



= 2V2t/c 



V2 = —=Aoa. 
V2 



(2.3) 



(2.4) 



(Note that all velocities in our model are taken to be 
"renormalized" velocities accounting for both band struc- 
ture and many-body effects within the context of Fermi 
liquid theory.) Utilizing these definitions, it becomes 
clear that the quasiparticle excitation spectrum in the 
vicinity of each of the gap nodes takes the form of an 
anisotropic Dirac cone 



Eh — 



where the degree of anisotropy is measured by the ratio of 
the two velocities. This ratio Vf /v2 appears prominently 
in the low temperature transport coefficients and is a 
measurable quantity which provides a convenient means 
of comparing theory to experiment. 

The low temperature physics of such a model is driven 
by the fact that at the four nodes, there is no gap to 
quasiparticle excitations. Hence, quasiparticles are gen- 
erated only in the vicinity of the gap nodes. This is very 
useful mathematically since it means that a momentum 
integral over the Brillouin zone can usually be replaced by 
a sum over nodes and an integral over the small region of 
k-space surrounding each node. Furt hermore, due to the 
form of the excitation spectrum (2.5), it is convenient to 
scale out the anisotropy of the Dirac cone, and change to 
polar coordinates in a new scaled momentum p = {p, 9) . 
Hence we will frequently make the substitution 



dkidk2 



E 



pdp 
2TrvfV2 



de 

2^ 



(2.6) 



A?. 



vjkl 



(2.5) 



where pi — Vfki — pcos9, p2 = = psin0, p = 

Vpi +P2 ^k, and po = ^/Wvjv^/a ~ 0(Ao) is a large 
scaled momentum cutoff defined such that the area of the 
new integration region is the same as that of the original 
Brillouin zone. Note that if quasiparticles are only gener- 
ated at the nodes and the rest of the Brillouin zone makes 
no contribution, then we should be safe in extending this 
limit to infinity. However, it is sometimes necessary to 
retain po through the intermediate stages of a calculation 
(usually as part of a ratio within a logarithm) maintain- 
ing throughout that all other energies are much smaller 
than this cutoff value. 

The fact that quasiparticles are concentrated in the 
vicinity of the gap nodes is also very useful when con- 
sidering impurity scattering. Since the initial and final 
momenta of a scattering event must always be approxi- 
mately equal to the k-space location of one of the four 
nodes, a general scattering potential, Vkk', need only be 
evaluated in three possible cases: intra-node scattering 
(k and k' at the same node), adjacent-node scattering (k 
and k' at adjacent nodes), and opposite- node scattering 
(k and k' at opposite nodes). These are depicted graph- 
ically in Fig. |l| and denoted respectively as Vi, V2, and 
V3. Hence, an arbitrary potential (varying slowly over 
the area of a node) is effectively reduced to a set of three 
parameters. This simplification proves quite helpful. 

Since the transport calculations in the sections that 
follow consist of the evaluation of Feynman diagrams via 
field-theoretic techniques, it is important to establish the 
types of Green's functions that will be utilized. For any 
superconductor, the existence of a condensate of ground 
state pairs means that the annihilation of an electron 
must be treated on the same footing as the creation of 
its mate, an electron with opposite momentum and spin. 
Hence we use the Nambu formalismEil in which the field 
operators are two-component spinors of the form 
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FIG. 1. Impurity scattering within d-wave model. Vi, V2, 
and V3 are the potentials for intra-node, adjacent-node, and 
opposite-node scattering. 



CfeT 



C-ki 



(2.7) 



and the resulting Green's functions are 2x2 matrices in 
Nambu space. Since we are concerned with finite temper- 
ature calculations, all diagrams will be evaluated using 
the Matsubara finite temperature formalisrr£3. Hence, 
the bare Matsubara Green's function expressed in Nambu 
formalism takes the form 



1 



El 



iu - efc 



(2.8) 



where the tilde denotes a Nambu space matrix and 
ioj = i(2n + is a fermionic Matsubara frequency. 

In the presence of impurities, the bare Green's function 
is dressed via scattering from the impurities and obtains 
a Matsubara self-energy, f^{iuj). Assuming that all but 
the scalar component of the self-energy can be neglected 
or absorbed into or A^, Dyson's equation yields that 
the dressed Matsubara Green's function is given by 



1 



Afe 



-El 



lUJ - 



Afc 
S(iw) 



efc 



(2.9) 



(Note that while this assumption has been explicitly jus- 
tified in both the Born and unitary scattering limits, the 
omitted self-energ^omponents can contribute for arbi- 
trary scattering. Ej'cfl For simplicity, we neglect such con- 
tributions in this investigation.) From the Matsubara 
functions, corresponding retarded functions are obtained 
by analytically continuing iuj ^ uj + i6 such that 



G'ret(k, w) = Gikjuj 



(2.10) 



and the impurity scattering rate is defined as 

r(w) = -lmS^et(w) (2.11) 

where 'Sreti^) — S(ici; lu + iS). 

With the Green's function in hand, it is a simple and 
illustrative step to calculate the density of states. In 
terms of the retarded Green's function, the density of 
states is given by 



N{uj) = - 



(2.12) 



Plugging in the Green's function ( |2.9| ), replacing the sum 
by a scaled integral about each node via (2^), neglect- 
ing the real part of the self-energy, and performing the 
integration we find that 



N{uj) 



1T^VfV2 



In 



Pa 

r(^) 




+ 



7TVfV2 



■ arctan 



r(c.)2 
2|w|r(w) 



(2.13) 



Note that in the absence of impurities {T{uj) = 0), 

N{uj) = (2.14) 

r(w)=o nvfV2 

while in the presence of impurities, there is a fpijte den- 
sity of quasiparticle states down to zero energyllj: 



iV(0) = 



-To In 



Po 



9 - i" 

TT-VfV2 1 



(2.15) 



where Fq = r(a; — > 0). These impurity-induced quasipar- 
ticles are responsible for the intriguing low temperature 
transport properties that we shall consider in the sections 
that follow. 



III. MICROWAVE ELECTRICAL 
CONDUCTIVITY 



Electrical conductivity can be calculated by means of 
the Kubo formulaE3 



a{n,T) = - 



(3.1) 



where Iiret{^) — n(zri — > VI + iS) and Ii{iQ) is 
the current-current correlation function (or polarization 
function) in the Matsubara finite temperature formalism: 



n(zf^)-- / dTe»"-(rjt(r)j(o)). 



(3.2) 



Thus, our first step is to derive an expression for the 
electrical current operator. Then by evaluating its corre- 
lation function, we obtain the electrical conductivity. 
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A. Electrical Current 

For a system of interacting electrons, the Hamiltonian 
is given by 



dxi/)t (x) 



2m* 



ijja (x) 



+ 1 [ dKdyi>iix)ijliy)Vi^~y)My)M^) (3-3) 



H 



where ipa (x) annihilates an electron of spin a at position 
X and V{x — y) is the electron-electron interaction po- 
tential. In the presence of a vector potential A(x), the 
Hamiltonian must be invariant under a local gauge trans- 
formation. While the second term is gauge invariant as 
written, the first must be modified by making the stan- 
dard replacement — iV — iV + eA. Thus, the Hamil- 
tonian becomes a functional of the vector potential and 
takes the form 



i/[A(x)] = / dxV4(x) 



2m* 



V'q(x) 



I I dxdy^t(x)V,t(y)y(x-y)V;0(y)V.„(x). (3.4) 



Note that only the kinetic term couples to the vector 
potential. Taking the functional derivative with respect 
to A(x) we obtain an expression for the electrical current: 



j^(x) - -Re 



5H 
(5A(x) 



2im 

Then taking the space-time Fourier transform yields 

e 

m 



(3.5) 



(3.6) 



and in the limit that q — > we obtain 

j^(0,l)) = -eEv/Vl/ti.I,,+^ 



(3.7) 



where v/ = dek/d\i ~ k/m* and we have expressed the 
final result in terms of 2 x 2 Nambu matrix notation. Note 
that in all cases, momentum indices on field operators 
denote both momentum and frequency. Although this 
result is well known, we have derived it here in order to 
provide a basis for comparison with the thermal current 
and spin current to be derived later. 



B. Bare Bubble (Electrical) 

Given the current, we proceed to calculate the corre- 
sponding current-current correlation function. The cor- 
relation function for the electrical current can be ex- 
pressed diagramatically as a fermionic bubble with fully 



dressed propagators and a fully dressed vertex in which 
each vertex contributes a coupling parameter e, a velocity 
V/, and a 2 X 2 Nambu formalism unit matrix 1. Assum- 
ing that the impurity scattering potential is isotropic in 
k-space, the corrections to the bare vertex vanish. Thus, 
in this approximation, the conductivity can be obtained 
from the calculation of a bubble with dressed propagators 
(i.e. Green's functions with self-energy included) but bare 
vertices (no interaction between the two propagators). 
Such a diagram will be referred to as a bare bubble. 

The calculation of the bare bubble polarization func- 
tion is of the same basic form for the electrical, thermal, 
and spin conductivities. Thus, to avoid repeating the 
same derivation several times, a generalized polarization 
function H^g"(ri) (applicable to all three cases) which 
depends on a coupling parameter g, a velocity v^, and 
a Nambu matrix Tq has been calculated in Appendix 
Applying the general result (A13) to the case of interest 
(5 = e, = V/, Tq = 1) we find that 



TT^ V2 J 27r 



duj 



00 

xTr G':,t{p,uj)G'Up,L^ + n) (3 



where 

Gret(p,t^) = 



{UJ - Erei(t^))^ -P^ 

P2 UJ - T,ret{^) - Pi 

In the universal limit (fi — > 0, T ^ 0), 
ni?(aj) — ni?(ti) + fi) dnp 



duj 



(3.9) 



(3.10) 



Hence, evaluating the rest of the integrand for a; ^ 0, 
noting that Sret(O) = — iFo, and integrating over momen- 
tum, we obtain the universal limit bare bubble electrical 
conductivity: 



(3.11) 



This is the universal conductivity obtained in Ref. |l} Fi- 
nite temperature corrections can be obtained via a Som- 
merfeld expansion (for r(ix')) and hojve been calcu- 
lated by Hirschfeld et alB and Graf et all 



C. Vertex Corrections (Electrical) 

The bare bubble conductivity derived above was calcu- 
lated in the approximation that vertex corrections could 
be safely neglected. It turns out (as we shall see later) 
that this approximation is justified if the impurity scat- 
tering potential is isotropic in k-space (Vkk' = V = 



4 



const). However, for a general scattering potential, cor- 
rections to the bare vertex can make a significant con- 
tribution and must be included in the calculation. To 
this end, we shall consider the contribution of the ladder 
corrections to the bare vertex (see Fig. ^ of Appendix 
. Once again, since the electrical calculation is of the 
same form as that for the thermal and spin conductivity, 
a generalized polarization function including vertex cor- 
rections (which can be applied to all three cases) has been 
calculated in Appendix ^ For the electrical conductiv- 
ity, the polarization function consists of a single bubble 
(with a dressed vertex) where each vertex contributes a 
coupling parameter e, a velocity v/, and a Nambu ma- 
trix 1. Plugging the se pa rameters into the generalized 
polarizat ion function ( B50| ) and using the electrical Kubo 
formula (3.1) we find that 



27r^ V2 J- 



duj 



n 



xRe 4°\uj,n)~ Jl'^'iuj^n) 



(0), 



(3.12) 



where Jf and are defined in Appendix In the 
universal limit (f7 ^ 0, T ^ 0), we can make use of 
(p^ to find that 



27r2 



V2 



Re 



r(0) 



'A2 ^2 



(0) 



(0) 

<A1 ^ lAl^l 



(3.13) 

where all functions are evaluated for Q,lu —^ 0. In these 
limits, the constituent functions defined in Appendix ^ 
take the form: 



F'(0) 



A'KVfV2 



i^Ji^ln^=4Ar(0) (3.14) 



7a 1 



V 



1 + (f iV(0))2X' 



1 + (f iV(0))2]/2 



ATrvfV2 



1T2 - {Al + Bl 



\n=l — \n=3 



2TTVfV2 



1^2 - 



(3.15) 

(3.16) 

(3.17a) 
(3.17b) 
(3.17c) 
(3.17d) 



/f^(0,0) = 



dF{z) 



dz 



= 2\n^-2 



/f(0,0)^I^^^21n^ 



/f^(0,0) 



^(0,0) = lim 



dF(z) 



dz 




Im [{uj — iTQ)F{LU + iTo 
2Fow 



(3.18) 



(3.19) 



-1 (3.20) 



1 (3.21) 



Thus, including vertex corrections, the universal limit 
electrical conductivity takes the form 



c^o 



TT^ V2 



(3.22) 



(0) 



'Bl 
'Al 



Inf 



27i°^nBi 



(0) 

Ta 



i 



(3.23) 



where Pvc is the scattering dependent vertex correction 
to the universal bare bubble result. Note that since the 
7's all depend on the difference between the intra-node 
T- matrix (n=l) and the opposite-node T-matrix (n=3), 
/3„c ^ 1 if the two scattering potentials are the same. 
Hence for an isotropic scattering potential, the bare bub- 
ble result (3.11) is recovered. However, in general we pre- 
sume that the scattering potential will fall off for large k 
and the potential for intra-node scattering will be larger 
than that for opposite-node scattering. If so, the 7's will 
be non-zero and the universal limit conductivity will de- 
viate from its bare bubble value. This correction to the 
conductivity due to differences between intra-node (for- 
ward) scattering and opposite-node (back) scattering is 
the node-discrete equivalent of the famous 1 — cos 9 factor 
obtained from vertex corrections in the conductivity cal- 
culation for a simple metal. As in the metallic case, the 
phenomenon at work is the fact that forward and back 
scattering can have different effects on the progress of a 
charge carrier. As a result, anisotropy in the scattering 
potential can renormalize the conductivity. 

In general, the evaluation of f3yc requires a numerical 
calculation since Fo must be obtained self-consistently 
as a function of impurity density and scattering poten- 
tial. Such calculations (presented in Appendix |^) indi- 
cate that for anisotropic scattering, the electrical vertex 
correction can be significant even to zeroth order in the 
impurity density. 

For the case of Born scattering, /3„c reduces to a more 
simple and illustrative form. In the Born limit (small V), 



5 



T„"(o) = \4i r,?(o) = 



(3.24) where /, 



^(0) _ (0) _ 



1^bI=0 (3.25) 



4TTVfV2 

and the zero-frequency scattering rate takes the form 



: po exp 



- 2Vi + Vi)N{Q) 



(3.26) 



where (as defined in Sec. |^ Vi, V2, and V3 correspond 
respectively to intra-node, adjacent-node, and opposite- 

— ^ 1 and defining 



node scattering. Noting that In ^ 



r,^27i%ln^.Jn.(T/,^ 
and a transport scattering rate 



Po 



Vi)N{0) (3.27) 



To-ri = -n,(2F22-f 21/32)iV(0) 



the vertex correction factor ( 3.23| ) reduces to 



i tr 



2Vi 



2V? 



2Vi 



(3.28) 



(3.29) 



kk' 



ifc, 



kk' 



The first term in brackets 



is peaked at = and is therefore a Fermi surface 
term (smeared over the extent of the gap). If we as- 
sume a circular Fermi surface, then replacing the k-sum 
by an integral in circular coordinates {k,6), presuming 
that Ak — A(0) <C Ep, and expanding the Landau func- 
tion in 2D harmonics 



1 



KO) 



J2 F! COS {(.{e - e')) (3.32) 



1=0 



this first term takes the form 



(3.33) 



where is the t — \ spin-symmetric Landau param- 
eter and J^(0) is the single spin normal state density of 
states at the Fermi surface. The second term in brack- 
ets is peaked bX — and is therefore a node term 
(contributing primarily at the gap nodes) . Replacing the 
k-sum by an integral about each of the nodes and noting 
by symmetry that 



4 

E 



OJ f 



(/ll /ll) 



(3.34) 



Note that the vertex correction depends on the scattering 
potential but is independent of the density of impurities. 
In this simple limit it is clear that if intra-node scattering 
is stronger than opposite- node scattering (as we expect), 
Tq will exceed Ttr and the universal limit electrical con- 
ductivity will be enhanced beyond the bare bubble result. 



D. Fermi Liquid Corrections (Electrical) 

To this point, our calculation of electrical conductivity 
has neglected the effects of the underl ying F ermi liquid 



we find that 



interaction between electrons. In Sec. Ill A , an explicit 
expression was derived for the electrical current in the 
absence of Fermi liquid interactions. In essence, it has 
the form 



Jo 



-6 2, V/fc Suko 



(3.30) 



where v/^ is the Fermi velocity at k and Suka is the 
deviation of the electron distribution from equilibrium. 
The Fermi liquid renormalization of such a current has 
been d erive d in Appendix |^. Plugging into the general 
result (Dll) we see that in the presence of Fermi liquid 
interactions, the electrical current is given by 



-ro 



e Suk' 



fc', 



k 



^fk fkk' 



A? 



2iel + Alf/^ 



El El 



^ 



arctan 



El 

Ek 



(3.31) 



h — Jo 



ro (/ll - /ll) 



47r2 



VfV2 



(3.35) 



where fl^ and f§i are respectively the intra-node and 
opposite-node spin-symmetric Fermi liquid interaction 
energies. In the small impurity density limit, the second 
(node) term can be neglected with respect to the first 
(Fermi surface) term. Hence, the renormalized electrical 
current is given by 



J =Jo"/i 
where (for a circular Fermi surface) 



fl 2 



To (/fl - /ll) 



Ff 
1 + f 



(3.36) 



(3.37) 



and the superscript 's' denotes that this is the spin- 
symmetric current renormalization factor. The simple 
form of this expression is due to our assumption of a cir- 
cular Fermi surface. For a more general Fermi surface, 
additional harmonics of the Landau function would have 
been generated. Thus, in practice, this factor should be 
treated as a parameter to be determined by experiment. 
Note that the renormalization applies to both the quasi- 
particle (normal) current and the supercurrent. However, 
since the supercurrent does not contribute to the real part 
of the AC conductivity, we are concerned only with the 
normal current. 

The basic physics of this renormalization is as follows. 
Upon the application of an electric field, quasiparticles 
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are perturbed to form a normal current. In the presence 
of the excited quasiparticles, the electron dispersion is 
modified due to the Fermi liquid interaction. The modi- 
fied dispersion yields a modified equilibrium distribution 
which means that the deviation from equilibrium is also 
modified and the current is renormalized. Note that the 
dominant term in the renormalization factor is a Fermi 
surface term resulting from the modification of the equi- 
librium condensate distribution in the presence of per- 
turbed quasiparticles. 

So far, we have discussed only how Fermi liquid interac- 
tions renormalize the electrical current density operator. 
Yet our goal is to determine the manner in which such 
interactions modify the electrical conductivity. In gen- 
eral, such modifications can be more complicated than 
merely renormalizing the constituent currents. However, 
as discussed in Appendix ^ current renormalization is 
the dominant effect in the T limit with which we 
are concerned. Therefore, since electrical conductivity is 
proportional to the current-current correlation function, 
two powers of our current renormalization factor appear 
in the conductivity. Hence, including both vertex correc- 
tions and Fermi liquid corrections, the electrical conduc- 
tivity in the universal limit takes the form 



<Jo = — Pvc afi' 



(3.38) 



Im 



7rsgn(w) [9{uj + po) ~ 9{uj - pq)] 
+ttpq[S{oj+po)-S{lj-po)]. (3.42) 



Performing the frequency integration yields that the nor- 
mal fluid density neglecting Fermi liquid corrections is 
given by 



/9"(r) _ 21n2w/ 



TT V2 



(3.43) 



which is precisely the result obtained in Rcf. ^ through 
an entirely different procedure. To include Fermi liquid 
corrections in the T — > limit, we note (via Appendix 
^ that the primary effect of Fermi liquid interactions 
is the renormalization of the current density operator. 
Thus, since the normal ffuid density is proportional to 
the current-current correlation function, we need only 
multiply by two factors of the current renormalization 
(obtained in the previous section) to find that 



p"(r) _ 21n2w/ 



TT V2 



(3.44) 



which agrees with the results of Refs. ^,|T^,^. Given this 
correspondence with prior work, we can be reassured of 
the accuracy of our calculations. 



IV. THERMAL CONDUCTIVITY 



E. Superfluid Density 

As a check on the accuracy of our conductivity calcu- 
lations, it is useful to make a brief digression and use our 
results to calculate an experimentally distinct quantity, 
the superfluid density (without impurities), p''(r). By 
definition (see Ref. @), 



p^(T)^p^(T = 0)-p"(T) 



(3.39) 



where p"(T) is the normal ffuid density. Hence, to obtain 
the temperature dependence of the superffuid density, it 
suffices to calculate the normal ffuid density. While the 
conductivity is related to the imaginary part of the polar- 
ization function, the normal ffuid density is proportional 
to the real part via 



p"(r) ReH,et(r! = 0) 



(3.40) 



Obtaining H^ef from the generalized result ( B49 ) in Ap- 
pendix setting r2 = 0, takin g the no impurities limit 
(r(tj) 0), and plugging into ( 3.40 ) we ffnd that 



p"(T) 1 vf 



m 



dcj npicj) Im l["\uj,0) (3.41) 



where 



Analogous to the case of electrical conductivity, ther- 
mal conductivity can be calculated by means of a thermal 
Kubo formulaE3: 



f 



1 ImH«,,(r!) 



rp2 



(4.1) 



where Tlreti^) = U''' (i^ ~* n + iS) and Il^iin) is the fi- 
nite temperature current-current correlation function (or 
polarization function). In this case, the appropriate cur- 
rent for the correlation function is the thermal current 
derived below. 



A. Thermal Current 

To derive an expression for the heat current in an 
anisotropic superconductor, we can follow the s-wave 
derivation of Ambegaokar and GriffinE a an d generalize 
to the case of an anisotropic gap. As in (3.3), the Hamil- 
tonian takes the form 



-V- 
2m* 



ijja (x) 



dxdyV^t (x)^Uy)F(x - y)V^^(y)Va(x) (4.2) 



Given the Hamiltonian, it is straightforward to obtain 
the equations of motion for the field operators 
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= ( 2^ + y dry (x - r)^t (j.)^^ (r) j (4.3) 

and to define a Hamiltonian density 
1 



2m' 



:VVl(x)-VVa(x) 



+ i / dyl^(x-y)V't(x)V'^(y)V/3(y)^a(x) (4.4) 



/i(x 



If all energies are measured with respect to the chemical 
potential, this Hamiltonian density is the heat density. 
Hence, the operator j'^(x) that satisfies the continuity 
equation 



/i(x) + V • jQ(x) = 



(4.5) 



can be interpreted as the heat current. Taking the time 
derivative of (4.4) and using the equations of motion ( |4.3| ) 
we find that 

/i = V • (^LVV'^.a + VV'LV'xa) dyV{y - x) 



(4.6) 



where the compact notation tl^xa = fpa (x) has been used 
f or t he sake of brevity. Defining j'^ = + we can use 
(15) to write 



(4.7) 



and 



V-j?(x)=.i J dyF(y-x) 



(4.8) 



Taking the space-time Fourier transform of ( [4.7[ ) in the 
limit that q ^ we obtain 



E 



(4.9) 



where v/ = dek/dk — k/m*, Cka is the space-time 
Fourier transform of '0a(x), and the second line is written 
in terms of the 2x2 Nambu matrix notation introduced 
in Sec. ||. Please note that in our compact notation, 
momentum indices on field operators always represent 
both momentum and frequency (i. e. '^k = ^(k, w) and 



Similarly, taking the space-time Fourier transform of 

(4.10) 



(4.8) generates four terms such that 

^q•j^(q,^7)=Xl+X2-yl-r2. 
The first such term, Xi , is given by 

k' s,uj' s 

X'5fe4-*:i-*:5-g'5fe3_fe2 + fc5(5i^^+cj2_cj3_[j^+Sl. (4.11) 

Taking the mean field approximation, retaining only the 
terms for which the average values are over (k t^^k I) 
pairs, and using the fact that (c[.|C^j.j^) is an even func- 
tion of uj , this becomes 



where 



'^k]'--ki 



(4.12) 



(4.13) 



Repeating this procedure for X2 , Yi , and Y2 and taking 
Afe to be real we find that 

q-j^(q,r!)=-^(Afc+g- Afe) 

X [w4^cl(^^_j^^ + {UJ + n)C-kiCk+q] (4.14) 
In the limit as q ^ 

d^k 



Ak + q - Ak 



q- 



9k 



q • V2 



(4.15) 



Thus, casting (4.14) in terms of the Nambu matrix for- 
malism we find that 



Yi 



V2*tf2* 



k+q 



(4.16) 



In the limit of small il, the second term can b e neglected 
compared to the first. Thus, combining (4.16) with (4.9) 
we obtain the following expression for the heat current 
in an anisotropic superconductor: 



vMfs-^k + q 



-V2^'j.fi1'fc 



(4.17) 



Note that for an s-wave superconductor, the gap is in- 
dependent of k and V2 = 0. T hus, the second term in 
the heat current vanishes and ( 4-L7 )-|reduces to the re- 
sult derived by Ambegaokar et alOcJ However, for a d- 
wave superconductor, the gap is anisotropic and V2 ^ 0. 
Hence, although the gap term may be small, neither term 
can be formally neglected. 
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B. Bare Bubble (Thermal) 



C. Vertex Corrections (Thermal) 



Unlike the electrical current, the thermal current has 
two terms: a "Fermi" term proportional to and and 
a "gap" term proportional to V2 and fi. Therefore, when 
we evaluate the current-current correlation function, we 
expect four bubbles rather than just one: Fermi-Fermi, 
Fermi-gap, gap-Fermi, and gap-gap. However, since the 
Fermi velocity and the gap velocity V2 are orthogo- 
nal at each of the gap nodes, the two cross terms cancel. 
Hence, the thermal conductivity has two terms: a Fermi 
term with velocity v / and Nambu matrix fs on each ver- 
tex and a gap term with velocity V2 and Nambu matrix 
fi on each vertex. For both terms the coupling parameter 
is {uj + n/2). 

Neglecting vertex corrections, each term can be ob- 
tained from the bare bubble generalized polarization 
function derived in Appendix Pluggi ng th e appro- 
priate parameters into the general result (A13) we find 
that 



k{VL,T) 




np{Lu) — np(LU + CI) 



In the universal limit (fi ^ 0, T ^ 0), 

nF{uj) — np{LU + Q) dnp 
Vl duj 



(4.18) 



(4.19) 



which for low T is very sharply peaked at lu — 0. Thus, 
evaluating the rest of the integrand for uj ^ Q, noting 
that Sret(O) — — iFo, performing the frequency integral 
via 



dnp 
duj 



duj 



(4.20) 



and integrating over momentum, we obtain the bare bub- 
ble thermal conductivity in the universal limit: 



T 



\ 1 vj + vl 

Kr 



VfV2 



(4.21) 



Neglecting the u| term in the numerator, this result and 
its finite temperature corrections, were originally calcu- 
lated by GratnCi alM The gap term was first obtained by 
Senthil et alB via a physical argument of Wiedemann- 
Franz correspondence with their expression for spin con- 
ductivity. It arises here as a direct result of the additional 
gap term found in our calculation of the thermal current 
for a d-wave superconductor. 



The bare bubble result derived in the previous section 
can be improved upon by including the contribution of 
the ladder corrections to the bare vertex (see Fig. ^ of Ap- 
pendix 1^) . A generalized polarization function including 
such vertex corrections has been derived in Appendix 
By pluggi ng th e appropriate parameters into this general 
formula (B50) both terms of the thermal conductivity 
(the Fermi term with parameters vj, fs, and lu + Q/2 
and the gap term with parameters V2, fi, and to + Q/2) 
can be obtained. Hence we find that 



K{n,T) 




°° np{u}) ~ np{LU + fl) 
duj 



2 r 



■ V2Re 



(4.22) 

where and are defined in appendix In the 
universal limit (fi — > 0, T ^ 0), the Fermi function factor 
is sharply peaked at w = 0. Thus, evaluating the J- 
functions for Q,lu 0, performing the frequency integral 
via ( 4.2Cl| ), and noting that 



(3), 



2 (0,0) = -jf)(0,0) 



J, 



jp^(o,o) = -j(^)(o,o) 



(4.23) 



we find that the universal limit thermal conductivity 
takes the form 



0^ = 

h'vc 



T 



1 v) + vl 

— / 

TT^ VfV2 



(4.24) 



(0) 
7a2 



(0) 

Tai 



y<0)(21n £2.- 



>A1 



21n£i)_2 



(4.25) 



where /3j^ is the thermal vertex correction factor and the 
7's are defined in ( ^.17 ). As for the electrical case, the 
thermal vertex correction must generally be evaluated 
numerically. The results of such numerical calculations 
(presented in Appendix |^) can be summarized as fol- 
lows: (1) For all scattering strengths (from Born to uni- 
tary) the thermal vertex correction is negligible compared 
to the electrical vertex correction. (2) In the small im- 
purity density limit, — 1 vanishes approximately as 



[ln(po/ro 



Thus, to zeroth order in the density of 



impurities, vertex corrections do not contribute. Hence, 



1 



(4.26) 
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and the universal limit thermal conductivity takes its 
bare bubble form 



1 vj+v^ 



VM^^^- (4.27) 

d / TT^ VfV2 

This is in stark contrast to the case of electrical conduc- 
tivity where we found a significant vertex correction even 
to zeroth order in the impurity density. 

D. Fermi Liquid Corrections (Thermal) 



As discussed in Appendix there may be additional 
corrections to the thermal conductivity due to underly- 
ing Fermi liquid interactions between electrons. In the 
T — > limit, the dominant effect of such interactions 
is to renormalize the current density operator. From 
(4.17) we see that the thermal current (in the absence of 
Fermi liquid interactions) has a rather complicated form 
including both a Fermi velocity term and a gap veloc- 
ity term. However, for the purposes of this analysis, it 
suffices to neglect the gap term (since it is known to be 
much smaller) and note that in essence, the bare thermal 
current has the form 



J? 



(4.28) 



V. SPIN CONDUCTIVITY 

For the spin conductivity case, the Kubo formula takes 
the form 



(5.1) 



where ILreti^) = W{in ^ n + iS) and is the fi- 

nite temperature current-current correlation function (or 
polarization function). Here the current that enters the 
correlation function is the spin current derived below. 



A. Spin Current 

To find an expression for the spin current operator in 
an anisotropic superconductor, we can write down the 
Hamiltonian and spin density operators and use the spin 
continuity equation to obtain the current. For spin, the 
continuity equation is 



p^(x) = -V.r(x) 



(5.2) 



where is the spin density and is the spin current 
density. The spin density equation of motion takes the 
standard form 



where Suka is the deviation of the electron distribution 
from equilibrium. To account for the effects of the Fermi 
liquid interactions, a general renormalization factor has 
been derived in Appe ndix ^ Plugging the appropriate 
parameters into (Dll) we find that the dressed thermal 
current is given by 



J ^ Jo 



k' a' 



Sn 



k' a' 



^k^fk fk. 



kk' 



El El + n 

~ Jkk' Jkk' 



-El 



arctan 



Ek 



(4.29) 



where ~ fH, + fH,. Since the entire summand is 
odd in Efc, the correction to the bare current integrates 
to zero. Hence, at least in the zero temperature limit, 
the thermal current is not renormalized by Fermi liquid 
effects. 



— Jo 



(4.30) 



Basically, due to the symmetry of the electron disper- 
sion about the Fermi surface, corrections to the thermal 
current cancel. Since the thermal current is unmodified, 
there are no Fermi liquid corrections to the universal limit 
thermal conductivity. Thus, including both vertex cor- 
rections and Fermi liquid corrections, the thermal con- 
ductivity retains the bare bubble form 



1 vj + vj 

TT^ VfV2 



(4.31) 



p^(x) = -z[p^(x),i7] 



(5.3) 



Thus, combining the two equations and using a Fourier 
representation for both the spin density and the spin cur- 
rent we find that 



(5.4) 



The Fourier transform of the spin density operator is 
given by 



k.u.a 
S 



(5.5) 



where = ±i and s = i. In the mean field approxima- 
tion, the Hamiltonian for a superconductor is expressed 
as 



H = 



-ki + c_fexCfe| 



(5.6) 



Thus, evaluating the commutator of (5.5) and (5.6) using 
fermionic anticommutation relations we obtain 



qj9== s 



i^k+q - Afe) [Af^-(k+q)l + C-kiCk+q] 



(5.7) 
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In the limit as 5 — s- 



A 



k+q 



Afc « q • 



q - V/ 



q - V2 



(5.8) 



(5.9) 



Hence, expressing the creation and annihilation operators 
in terms of 2x2 Nambu matrix notation we find that 



k.uj 



(5.10) 



Note that the spin curr ent t akes precisely the same form 
as the thermal current (4-17) with an appropriate change 
of coupling parameter. 



B. Bare Bubble (Spin) 

As in the thermal conductivity case, the spin current 
has both a "Fermi" term and a "gap" term. Hence, evalu- 
ating the current-current correlation function and noting 
(as before) that the cross-terms cancel, we find that the 
spin conductivity consists of two bubbles: a Fermi term 
with V/ and fa on each vertex and a gap term with V2 
and fi on each vertex. These are precisely the bubbles 
that we evaluated for the thermal case except that here 
the coupling constant is the spin (s — 1/2) rather than 
the frequency. 

Neglecting vertex corrections, each of the two bubbles 
can be evaluated by plugging the appropriate set of pa- 
rameters into the bare bubble generalized polarization 
fimction (A13) derived in Appendix Doing so we find 
that 



■n^VfV2 



271 



dco 



np{Lu) — np{uj + f2) 



vjTi 



^ ret 



g;'^,(p, w)fiG;'^,(p, uj + ^)fi 



In the universal limit (il ^ 0, T — > 0), 

npioj) — npi'^ + ^) dnp 



duj 



(5.11) 



(5.12) 



Thus, evaluating the rest of the integrand in the w — *■ 
limit, noting that l]ret(0) = — jFo, and integrating over 
momentum, we obtain the universal limit bare bubble 
spin conductivity: 



s2 + vl 
TT^ VfV2 



(5.13) 



This agrees (aside from a disputed factor of 2) with the 
result obtained by Senthil et alB 



C. Vertex Corrections (Spin) 

As in the electrical and thermal cases discussed previ- 
ously, the bare bubble result derived above can be im- 
proved upon by including the contribution of the ladder 
corrections to the bare vertex. By plugging the appro- 
priate parameters into the generalized polarization func- 
tion (including vertex corrections) derived in Appendix 
both the Fermi term and the gap term of the spin 
conductivity can be obtained. Doing so wc find that 



2lT^VfV2 



duj 



np{ijj) — np(uj + Vl) 



where Jf and are defined in Appendix B3 
unive rsal limit (il 0, T ^ 0), we can ma 
( ^.12|) to evaluate the frequency integral and find that 



(5.14) 

In the 
ce use of 



VfV2 



-0" 



(5.15) 



where the spin vertex correction factor, is identical 
t o th e thermal vertex correction factor, defined in 
(4.25). Thus, mirroring the analysis described in Sec. 
[VC, we note that for small impurity density, (3^^ ~ 1 
and the universal limit spin conductivity takes its bare 
bubble form 



VfV2 



(5.16) 



As in the thermal case and in contrast to the electrical 
case, vertex corrections do not contribute (to zeroth order 
in the impurity density). 



D. Fermi Liquid Corrections (Spin) 

An additional correction to the spin conductivity is 
needed to account for the effects of underlying Fermi liq- 
uid interactions between electrons. In the T — + limit, 
the dominant effect of Fermi liquid interactions is the 
renormalization of the current density operator (see Ap- 
pendix ^). Neglecting the gap term in (5.1C) (since it is 
small and rather difficult to deal with), the spin current 
has the basic form 



Jo = / . SaVfk Sn 



ka 



(5.17) 



where Sa = ±1/2. The Fermi liquid renormalization of 
this current can then be obtained by pluggi ng th e appro- 
priate parameters into the general result (Dll) derived 
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in Appendix Noting that the math is completely anal- 
ogous to that of the electrical calculation in Sec. [II D 



we can easily adapt the electrical result to the present 
case. Replacing — e with s (= 1/2) and changing spin- 
symmetric designations to spin-antisymmetric ones, we 
find that the Fermi liquid renormalization of the spin 
current takes the form 



(5.18) 



ization factor which, for a general Fermi surface, is some 
complicated function of the spin-antisymmetric Landau 
parameters, F^. For the simplified case of a circular 
Fermi surface. 



a 



fi 



1 



pa 



(5.19) 



As in the electrical case, the current renormalization 
is dominated by a Fermi surface term resulting from the 
interaction-induced modification of the equilibrium dis- 
tribution of the condensate. At first glance, this result is 
a bit surprising since the Zeeman field which generates 
the (normal) spin current cannot induce a supercurrent. 
However, it must be understood that the renormaliza- 
tion of the normal current has nothing to do with the 
presence of a supercurrent. Rather, due to the existence 
of the normal current, the equilibrium distribution of the 
condensate is modified (via Fermi liquid effects). It is 
this modification of the equilibrium condensate, not the 
presence of an excited condensate (supercurrent), that 
gives rise to the renormalization of the spin current. 

Since the spin conductivity is proportional to the 
current-current correlation function, it obtains two fac- 
tors of the spin current renormalization. Thus, including 
both vertex corrections and Fermi liquid corrections, the 
universal limit spin conductivity takes the form 



VfV2 



(5.20) 



VI. CONCLUSIONS 

In the presence of impurities, the gap symmetry of a 
d-wave superconductor yields the generation of impurity- 
induced quasiparticles at the gap nodes. The transport 
properties of the resulting system are quite unique since 
such quasiparticles are both generated and scattered by 
impurities. In the — > 0, T ^ limit, bare bubble 
calculations indicate that transport coefficients are "uni- 
versal" , independent of the impurity density or scattering 
rate. However, once the contributions of vertex correc- 
tions and Fermi liquid corrections are included, we find 
that (putting in the ?i's) the electrical, thermal, and spin 
conductivities in this universal limit take the form 



Co 

Kq 

T 



s 



'Hi 

V2 



V2_ 



V2_ 

Vf 



(6.1a) 

(6.1b) 
(6.1c) 



where a?; is the spin-antisymmetric current renormal- (E^) ^^'^ 



where P^c is a scattering dependent vertex correction 



and a'j:i 



are spin-symmetric and spin- 



antisymmetric Fermi liquid factors ( ^.36| . 5T8| ). Note that 
these are the 2D conductivities of a single Cu02 plane. 
To obtain 3D conductivities, they must be multiplied by 
the number of Cu02 planes per unit length stacked along 
the c-axis. 

The "law" of Wiedemann and Franz suggests that the 
transport coefficients should be related such that 



K 

of 



3 e2 



K 



3 s2 



(T 

a 



(6.2) 



However, examination of the expressions above yields 
three sources of Wiedemann-Franz violation: current 
operator definition corrections, vertex corrections, and 
Fermi liquid corrections. First of all, since the electri- 
cal current has only a Fermi term while the thermal and 
spin currents include both a Fermi term and a gap term, 
(To is proportional to the ratio Vf /v2 while ko and CTq in- 
volve an extra V2/vf term. These extra terms arise when 
the thermal and spin current operators are corrected to 
account for the anisotropy of thej-arder parameter. How- 
ever, since Vf /v2 ~ 14 for YBC0t2l, this type of violation 
is of more qualitative than quantitative importance. Sec- 
ondly, unless impurity scattering is completely isotropic 
in k-space, the electrical conductivity contains a scat- 
tering dependent vertex correction, f3vc, which cannot 
be neglected even to zeroth order in impurity density. 
However, analogous corrections to the thermal and spin 
conductivities vanish in the small impurity density limit. 
Thus we expect a scattering dependent enhancement of 
CTo that is absent in kq and CTq. Finally, due to underlying 
Fermi liquid interactions, the electrical and spin conduc- 
tivities gain spin-symmetric and spin-antisymmetric cor- 
rection factors respectively. Corresponding corrections 
to the thermal current cancel due to particle-hole sym- 
metry. Hence, while Fermi-liquid interactions modify co 
and ctq, the value of kq is unaffected. 

The physical origin of the first two corrections lies 
with the velocity dependence of the current operators. 
Although somewhat obscured in the Nambu formalism, 
when our current operators (3.7,4.17 5.1C) are rewritten 
in the quasiparticle basis, it is clear that the electri- 
cal current is proportional to the Fermi velocity, Vf — 
dek/dk, while the thermal and spin currents are propor- 
tional to the group velocity, vg = dEk/dk. This differ- 
ence arises because quasiparticles carry definite energy 
and spin but do not carry definite charge. Since energy 
and spin are well defined in the quasiparticle basis, ther- 
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mal and spin currents are proportional to the group ve- 
locity, the derivative of the quasiparticle dispersion. By 
contrast, the electron and hole parts of each quasiparti- 
cle have opposite charge and opposite velocity. Therefore 
each part carries the same electrical current, proportional 
to the normal state Fermi velocity. This point was em- 
phasized in Ref. For a d-wave superconductor where 
both efc and are momentum-dependent, the group ve- 
locity will have both a v/ component and a V2 component 
while the Fermi velocity can only have a v/ component 
(see Fig. ^). This is the source of the extra gap terms in 
the thermal and spin conductivities. (Similar conclusions 
were drawn by Moreno and Coleman.tj) 

The role of vertex corrections can be understood by 
considering the graphical depictions of the Fermi veloc- 
ity and group velocity presented in Fig. |^. Throughout 
the area of a node, the magnitude and direction of the 
Fermi velocity is approximately constant. Thus, the elec- 
trical current can relax much more effectively via scatter- 
ing from node to node than it can via scattering within 
a single node. It is therefore necessary to distinguish, 
mathematically, between the effects of intra-node scat- 
tering and inter-node scattering. This is accomplished 
through the inclusion of vertex corrections. In contrast, 
the group velocity varies significantly over the area of 
a node. Therefore, the thermal and spin currents can 
relax through cither intra-node scattering or scattering 
between nodes. As a result, the different types of scat- 
tering play nearly the same role and need not be distin- 
guished. Hence, vertex corrections do not contribute to 
the thermal and spin conductivity. 

Fermi liquid corrections result from the redistribution 
of equilibrium electrons in response to the presence of 
interactions between excited electrons. In essence, this 
redistribution gives rise to a drag current that can renor- 
malize the quasiparticle current and therefore the con- 
ductivity. The character of the renormalization depends 
on the nature of the coupling parameter for a particu- 
lar current. Since the spin current gets opposite contri- 
butions from the two species of spin, the spin conduc- 
tivity gets a spin-antisymmetric renormalization. How- 
ever, charge is spin-independent so the electrical conduc- 
tivity gets a spin-symmetric renormalization. Further- 
more, since energy changes sign across the Fermi surface, 
particle-hole symmetry dictates that the effects of Fermi 
liquid interactions on the thermal current must cancel. 
Thus, thermal conductivity is not renormalized. 

The velocity ratio, Vf/v2, is a fundamental material pa- 
rameter which measures the anisotropy of the quasiparti- 
cle excitation spectrum. Therefore, an important objec- 
tive in measuring quantities such as the normal fluid den- 
sity and the universal limit transport coefficients, which 
all depend on Vf/v2, is to obtain the value of this ra- 
tio. However, due to vertex corrections and Fermi liquid 
corrections, the electrical conductivity, spin conductivity, 
and normal fluid density depend on parameters (such as 
interaction energy and/or scattering potential) with val- 
ues that are not well known. Only the thermal conduc- 







(b) 







FIG. 2. Schematic depictions of the (a) electrical current 
and (b) thermal/spin current in the vicinity of the four gap 
nodes. Electrical current is proportional to the Fermi velocity, 
V/ = 9efc/9k, whereas thermal/spin current is proportional 
to the group velocity, vg = dEk/d\i.. The ellipses drawn at 
each node denote (at a very exaggerated scale) the regions of 
k-space within which impurity-induced quasiparticles are gen- 
erated in the universal limit. (Note that in the small impurity 
density regime with which we are concerned (Fo -C Ao), these 
nodal regions are point-like on the scale of the Brillouin zone.) 



tivity involves neither vertex corrections nor Fermi liquid 
corrections. Thus, kq is the only truly "universal" coef- 
ficient and is the quantity from which the value of Vf /v2 
can be most directly obtained. On the other hand, the 



linear T coefficient of the superfluid density (3.44) is pro- 
portional to ^ Vf /v2- Hence, these two measurements 
can be combined to determine the Fermi liquid factor, 

In jfact, while this paper was in preparation, Chiao 
et alx3 applied these conclusions to the results of a se- 
ries of recent experiments performed on optimally-doped 
Bi2Sr2CaCu208 (BSCCO). By analyzing the residual 
linear term in their very low temperature thermal con- 
ductivity measurements in terms of (6.1b), they obtained 
a value for the velocity ratio, Vf/v2 — 20. This is the 
same value obtained from the ARPES measurements of 
Mesot et aLlla Going further, by combining this result 
with the linear T coefficient of the su«erfiuid density mea- 
sured by Waldram and co-workerscZI and making use of 
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44), they extracted a value for the Fermi-hquid correc- 
tion, ^ = 0.41. These results provide an experimental 
verification of our analysis. 
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APPENDIX A: BARE BUBBLE CALCULATION 

In the absence of vertex corrections, our calculations 
of electrical, thermal, and spin conductivity all require 
the evaluation of bare bubble diagrams depicting various 
types of polarization functions. The details of these dif- 
ferent calculations are all quite similar. They differ only 
in the coupling parameter, velocity, and Pauli matrix 
contributed by each bare vertex. Rather than repeat- 
ing the same basic calculation several times, it is conve- 
nient to calculate a generalized polarization function here 
which can be referred to for each of the specific cases of 
interest. This generalized function H^^", will depend on 
a coupling parameter g, a velocity v^, and a Nambu space 
Pauli matrix fa where 

g = {e,s(= l/2),w + f7/2} 
■V( = {v/, V2} 
fa = {rai^ i),fi,f2,f3}. 

Evaluating the diagram in Fig. ^ we find that 



n9^"(^f]) = -Y^g^ve^reTT\g{k,^Lu)fagik,iu; + in)fa 



(Al) 



where Q (k, iuj) is the 2x2 Nambu matrix form of the Mat- 
subara Green's function. For a d-wave superconductor at 
temperatures much less than the gap maximum, quaa: 
particles are generated primarily at the four gap nodes Jj' 
Thus, linearizing the quasiparticle spectrum about the 




nodes and defining a coordinate system (^1,^2) at each 
node with ki (k2) perpendicular (parallel) to the Fermi 
surface, we can replace our momentum sum by an inte- 
gral over the k-space area surrounding each node. If we 
further define a scaled momentum {pi,P2) we can let 



k j=i 



_fk_ 

(2^ 



4 

E 



{2TT)^VfV2 



(A2) 



where pi = Vfki = and p2 = ^2^2 = ^k- Since 
V/ = Vfki and V2 — W2k2 at each of the four nodes, the 
sum over nodes yields 



(A3) 



This, in turn, allows the definition of a scalar polarization 
function via 

n5^"(if7) = nf^"(if^)a. (A4) 

Defining a spectral representation for Q, we can write 



-duJi 



where 



A{p,uj) = --G;',,(p,c^) 

TT 



(A5) 



(A6) 



and G'/gf is the imaginary pa rt o f the retarded Green's 
function. Plugging back into (Al) we obtain 



A{p,UJi)faA{p,UJ2)fa 



xTr 



S 



where 



1 



1 



9 



uji iuj + ifl ~ UJ2 



(A7) 



(A8) 



Evaluating the Matsubara sum in the standard wayO we 
pick up a contribution from each of the poles of the sum- 
mand. Since the intermediate results differ depending on 
the frequency-dependence of the coupling parameter, it 
is best to handle the frequency-independent coupling and 
frequency-dependent coupling cases separately. 

For g = {e, s} (frequency- independent coupling), the 
sum is straightforward. Adding the contribution of the 
two poles and then continuing ifi ^ + «(5 we obtain the 
retarded function 



FIG. 3. Bare polarization bubble. 



Sret 

where 



ciJi — C1J2 + " + io 
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np{ijj) 



(AlO) 



is the Fermi function. 

For g = ll) + (frequency-dependent coupling), we 
proceed in the same way but there are a few technicahties 
that must be clarified. First of all, it should be under- 
stood that within the Matsubara sum we really mean 
g — > iuj + in/2. Only after the sum has been evalu- 
ated and all frequencies have been continued to the real 
axis should the stated form of the coupling parameter be 
taken literally. Secondly, note that with this frequency 
dependent g, the summand has two extra powers of fre- 
quency. As a result, the sum appears to be divergent. 
However, as discussed by Ambegaokar and GriffinE3, this 
apparent divergence results from an improper treatment 
of time-derivatives within the time-ordered correlation 
function and should be ignored. Doing so, we proceed 
just as before. Adding the contribution of the two sum- 
mand poles and continuing the external frequency to the 
real axis we find that 



LOi — UJ2 + n + i5 



(All) 



Plugging (A£) and ( All ) back into (A7), writing the 
spectral function in terms of the retarded Green's func- 
tion via (A6), and using the identity 



1 1 

— P moix) 



(A12) 



to take the imaginary part, we find that 
1 vj f (Pp 



Imn?^"^ow 



TT^ VfV2 .1 2n 
2 



duj 



xg {npi^^ + n) — npiuj)) 
xTT[G'Up,Lo)fo.G':,tip,u + n)fo] (A13) 

for all three coupling parameters g = {e, s,uj + il/2}. Ne- 
glecting vertex corrections, this is the imaginary part of 
the generalized retarded polarization function. The real 
part can be obtained via Kramers-Kr onig analysis. By 
specifying different input parameters, (A13) can be used 
to obtain the electrical, thermal, and spin conductivity. 



APPENDIX B: VERTEX CORRECTIONS 




(C) 



+ 



+ 



\ + 



FIG. 4. (a) Polarization bubble with dressed vertex, (b) 
Ladder series, (c) T-matrix series. 



that = vi'k.i we find that the generalized polarization 
function takes the form 



iu k 



xTr 



^(k, iw)f„^(k, iuj + irj)faf ^"(k, iu, in) (Bl) 



which is equivalent to the bare bubble result ( |Al| ) with 
the unit vector k^ from the second bare vertex replaced 
by a more general vertex function F^". Evaluating the 
diagram series in Fig. §(b) we obtain an equation which 
can be solved for the vertex function: 

f„f ^"(k) = kefo, + n,Y, Tkk" {iuJ + in)g{k", iu; + iil) 

k" 

X tJ:'" (k")^ (k" , iuj)fknk{ioj) (B2) 

where is the impurity density and Tkk' {i^) is the im- 
purity scattering T-matrix defined by the diagram series 
in Fig. 0(c). Multiplying from the left by Tq, we can define 



f^"{k,iuj,in)=\<ie l+k'iiu^in) 



(B3) 



where 



k^A" = Tii f^Tkk" {i^ + 'i^)G(k" , iu + iQ)fa 



Unless the scattering potential is completely isotropic, 
the bare bubble results of Appendix ^ can be improved 
upon by including the contributions of vertex corrections. 
In this section the ladder corrections depicted in Fig. ^ 
will be included. Once again, our object is to obtain an 
expression for a generalized polarization function, 11^^", 
in which each of the vertices contribute a coupling pa- 
rameter g, a velocity v^, and a Nambu space Pauli ma- 
trix fa- Evaluating the diagram in Fig. ||(a) and noting 



xT"^{k")g{k",iu;)Tk"k{tu;). (B4) 



Acting on both sides of ( |B3[ ), using (B4) to replace the 
left hand side by k^A", and noting that by the symmetry 
of the scattering potential 



^k,^k^^(k'.k) 



(B5) 



we find that 
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A" = ^ (k' • k)Tafk'k{i(^ + in)^(k, icu + iQ)Ta where 



X a + A" 



g{'k,iuj)Tkk'iiuj) 



(B6) 



where k' is the final momentum of a scattering event. 

An expression for the T-matrix in terms of the scatter- 
ing potential Vkk' and the Matsubara self-energy S](zcj) 
(itself a function of the scattering potential) can be ob- 
tained by evaluating the diagram series in figure ^(c). 

fkk'{ioj) = hVkk' + ^T3^(ki,ia;)f3T4fciVfeife' H 

ki 

(B7) 

As discussed in Appendix at temperatures low com- 
pared to the gap maximum, quasiparticles are only gen- 
erated in the small regions of k-space surrounding each 
of the four gap nodes. Hence we can let 



cPp 



(B8) 



and note that the initial and final momenta of a scat- 
tering event must always be approximately equal to the 
k-space location of one of the four nodes. Thus if j and 
j' are node indices (1 to 4), the scattering potential takes 
the form 



kk' 



(B9) 



where T£ is a 4x4 matrix in node space. Due to the sym- 
metry of the nodes, V_ consists of only three independent 
parameters: Vi for intra-node scattering, V2 for adjacent- 
node scattering^ and V3 for opposite-node scattering, (see 
Fig. in Sec. 0) 



V 



(Vi V2 P3 V2 

Vi Vi V2 

V2 Vi V2 

\V2 V^3 V2 Vx 



(BIO) 



Hence ( |B7[ ) becomes 



-T3 



(2n)^VfV2 



(Bll) 



Performing the node integration we see that the inte- 
gral of the Green's function about a node is a scalar in 
Nambu space (proportional to 1). Since raised to an 
even power is equal to 1, the T-matrix splits into a f3 
component and a 1 component. Summing the resulting 
geometric series we find that 



33 



V 



jj' 



rph 

'jj' ^ 



1 - i?'(zc^)2y V „v 



(B12) 



F'{iuj) = 



F{iuj - S(iw)) 



'i'KVfV2 



F{x) = xln 1 



(B13) 



(B14) 



and po is the large scaled momentum cutoff defined in 
Sec.|l| 



Taking the Nambu space trace of (BE), cyclically per- 
muting within the trace, and replacing the mom entum 
sum with a scaled integral about the nodes via (p3^) yields 



TrA" = Trk / —^^g{vM 

4 



;g(p,iw + ir2)f„ (^i + A") . (B15) 



Defining node j' to be node 1 we can write for j = 
1,2,3,4 



kJ' -k^ ^ {1,0,-1,0} 

{fi,f2,f3,f2} 



^33 ~ ^3 3 



and therefore 
4 



^k^ • k^J Tjji {iLu)faTj'j{iLu + iQ) 



'i7TVfV2 

m 



T3 T, 



(B16) 



where 



lA 



iB 



'i'KVfV2 



n=3 



4:TTVfV2 



1 - |n=3 



and we define 



+1 for a = 0, 3 
-1 for a = 1,2 



(B17) 



(B18) 



(B19) 



Plugging ( B16 ) into ( B15 ) and defining 



I°'{iu3,iVL) 



Q{-pTiu})TaQ{-p,iu) + in)Ta (B20) 
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and 

we find that 
Tr A" = Tr 



g{p,iuj)f3fag{p,iuj + iQ)fa (B21) 



72/"(i + A") 



Tr 



7g/3"(i + A") . 

(B22) 



Similarly, multiplying (B6) by and repeating our steps 
we find tliat 



Tr 



T3A^ 



= Tr 



7g/"(i + A") 



Tr 



7S/3"(i + A") . 

(B23) 



At this point it is useful to carry out the momentum 
integrals in (B20|JB2T]) to obtain an explicit form for /" 
and Recalling the form of the Matsubara Green's 
function from Sec. ||, noting that pi ~ et and p2 = A^, 
and defining 



we can write 



g{p,iu;) = 



fit +P1T3 +P2T1 



(B24a) 
(B24b) 



(B25) 



Further, using the definition of r]^ from (B19) and simi- 
larly defining 



-1 for a = 0, 1 
-1 for a = 2,3 



(B26) 



we find that 



Tay[p,lUJ + lil)Ta ^ -5 5 • (B27) 

J2 ^ P 



Thus, plugging ( |B25D and (|B27| ) into ( |B20D , noting that 
Pi = pcosO and p2 — psm9, and performing the angular 
integral we find that 



J" = = i r 2py^/i4i^^^dp (B28) 
Jo 



' if? - p')ifi - p') 



where 



+ 1 for a = 

for a = 1 , 3 
-1 for a = 2 



(B29) 



Factoring the integrand and performing the p-integral 
yields 



(/i + a^f2)F{f2) - (/2 + a^h)F{h) 



(B30) 



where F{x) is defined via (B14). Similarly, acting on 
(p2l|) yields that 



Jo 



if!-p')ifi^p') 



where 



+ 1 for a = 3 

for a = 0, 2 
— 1 for a = 1 



(B31) 



(B32) 



It is easy to see that a'^ = aa+i (where the index addition 
is defined modulo 4). Hence 



(B33a) 
(B33b) 



where /" is given by (B30). 

Now that /" and /" have been evaluated, ( B22| ) and 
( B23 ) become a set of coupled equations for TrA" and 
Tr[f3A°]: 

Tr A" = 73/"(2 + Tr A") + 7g/"+iTr[f3A"] (B34a) 
Tr[f3A"] ^ 7g/"(2 + Tr A") + 72/"+^Tr[f3A"] (B34b) 

Solving simultaneously yields 

2 



Tr 



1 + A" 



1 - 



73^" ( 



1 



7g 7g^° + ' 
74 1-7^/° + ^ 



(B35) 



This is a v ery us eful result since using (B3), (| 
and (^H) with (^) yields that 



(B20) 



i + A" 



(B36) 



Thus plugging (B35|) into ( |B36D and making use of ( |A3| ) 
from Appendix A| we find that 



where 



J" 



vf 1 



nvfV2 13 



I" 



1% + ' 



1-7^/" l + ^T^ 
'A \ 7° 1-7°/ 



(B37) 



(B38) 



(B39) 



and we note that 7"+^, 7^, and 7^ are all functions 
of iuj and iU,. 

Provided the input self-energy is of a proper func- 
tional form, J"(z,ir2) will be analytic throughout the 
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complex plane except for two branch cuts at Im z 
and Imz = — 17. Thus, evaluating the Matsubara suniE3 
we pick up a contribution from each of the branch cuts 
of the summand. Consequently, it is useful to consider 
the form of J"(2;, ifi) above and below each of the branch 
cuts. Upon exami nation of the frequency dependence of 
this function via ( p339| ,B17-B21) it is clear that the in- 
ternal and external frequencies, iio and iil, enter only 
through functional couplets of the form 



P(icj, i^) = A(iLS)B{iijj + i^). 



(B40) 



Furthermore, due to the defined analytic structure of 
Matsubara Green's functions, the functions composing 
these couplets always have a Matsubara-like analytic 
structure and satisfy 



(B41) 
(B42) 



Consider the form of such a couplet above and below the 
branch cuts of our summand. Defining 



Pi{ijj,in) 



P{uj 
P{uj 
P{lo 
P{lo 



iS, ifl) 
iS, ifl) 
iQ + iS, ifl) 
ifl — iS, ifl) 



and continuing ifl fl + iS we see that 

P3{u;,n) = P2{Lj-n,fi) 



(B43a) 
(B43b) 
(B43c) 
(B43d) 

(B44a) 
(B44b) 



Since J" is composed of such couplets, if we similarly 
define 



j^{uj,in) EE .riLj-iS,^) 

it follows that 



(B45a) 
(B45b) 
(B45c) 
(B45d) 

(B46a) 
(B46b) 







These relations will be very helpful in what follows. Now 
we can proceed with the Matsubara sum. As in Appendix 
it is best to treat the frequency-independent coupling 
and frequency-dependent coupling cases separately. 

For g — {e,s} (frequency-independent coupling), the 
sum is straightforward. Adding the contribu tions of the 
two branch cuts and using the definitions in ( B45 ) yields 



9 n 



27r2 VfV2 J_r 



(hjnp{ijj) 



(B47) 



where np^Lo) is the Fermi function. 

For g = u> + fl/2 (frequency-dependent coupling), not- 
ing the technical issues discussed in the analogous stage 
of the bare bubble calculation (see Appendix |A|) , we can 
proceed as above. Ad ding t he contributions of the two 
branch cuts and using ( B45 ) we find that 



1 vj 
27r^ VfV2 



duj np{uj) 



U) 

2 



{J^{uj,iVL)- J1{uj,iVC)) . (B48) 



In either case, continuing zil ^ 17 -I- «(5, making use of 
the relations in (B46), and shifting w — > w -|- 17 in the last 
two terms, we obtain the retarded polarization function 



1 vi 



duj g^ 



27r2 VfV2 
nF{uj){J'^{uj,n) - J2"(w,17)) 

np{u; + n){jfiuj,n) - J^{uj,n))] (B49) 



which is valid for all three coupling parameters. Tak- 
ing the imaginary part, noting that Re[z] = Re[z*], and 
expanding Jf and J2 yields 



Imnf,^(17) 



xRc 



27r2 VfV2 J- 



duj g^ (np{LO + ri) — np(uo)) 



-'2 



Ta2 1-722-^2 



I? 



(B50) 



where 



jc, _ (/l+«a/2)F(/2) - (/2 + «c/l)f (/l) 

J2 Jl 



.0. o^ _ (/r+aa/2)f (/2) - (/2+aa/r)f (/r) 

J2 ~ Jl 

fl=UJ- S,-et(t^) 
/2 = o; + 17 — Yireti^ + 



(B53a) 
(B53b) 



lAl 



AllVfV2 



+r„^(c^)r„^(c. + 17)) [|„=i 



|n=3 



(B54) 
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7A2 



4:TrvfV2 



7s 1 



1 — |n=3 



+t:{uj)t::{co + n)) [I 



n— 1 |n— 3 



(B55) 



(B56) 



to make all quantities dimensionless by expressing en- 
ergies in units of po and lengths in units of 2a. This 
choice of units sets the frequently encountered constant, 
A'KVfV2, equal to 1. For a particular set of input pa- 
rameters, the computation is done in two steps. First, 
solve self-consistently for the zero-frequency scattering 
rate, Fq, via 



(C3) 



1B2 



A7rvfV2 



1 - |n=3 



V 



-F'{lo)V^ 



F'iu;) = 



F[UJ - Y,ret[uj)) 
ATTVfV2 



(B57) 



(B58) 



(B59) 



(B60) 



and F{x) is defined in (B14). The above equations define 
the imaginary part of the generalized retarded polariza- 
tion function including ladder corrections to the vertex. 
By specifying different input parameters we can use it to 
obtain the electrical, thermal, and spin conductivity. 



APPENDIX C: NUMERICAL ANALYSIS OF 
UNIVERSAL LIMIT VERTEX CORRECTIONS 



In Sees. [lie, IV C, and VC, the vertex correction 



factors for the universal limit electrical and thermal/spin 
conductivities were found to be 



which (in the universal limit with our choice of units) 
reduces to 



2ri, In ■ 



i^l+(2rolni3-) V')^^ 



1. (C4) 



Then given Fq, plug into (CI) and (C2) to obtain the 
vertex corrections. 

The result will depend on the set of four input parame- 
ters which determine the impurity density and scattering 
potential, {n^, Vi, V2, V3}. In our units, rii = 4z where z 
is the substitutional fraction of impurities. Furthermore, 
it is convenient to parametrize the scattering potential 
via 



Vi = Vs, 



. tan 



R2 



V2 



V3 

^ (C5) 



where Vgcaie is a dimensionless constant and d ranges 
from to 1. Hence, the scattering anisotropy is given by 
i?2 and i?3 while the scattering strength is given by d. 
Note that while is tempting to think of as some sort 
of scattering phase shift, to do so would be stretching 
an analogy beyond its realm of usefulness. The above 
is merely a helpful parametrization, the aim of which is 
to allow us to go smoothly from the Born limit to the 
unitary limit as d varies from to 1. For our purposes, 
these limits are determined by the denominator in {C4): 



Pvc — 



w;ln^(ln^-l) 



l-27i°^n^ 



1 - 2 



^(0) _ -r'ni Win £0 _ 1 



(CI) 



1 - -7^°' 



l + 7i"^ 



1 ^^7!{\'l-7i">(21n^_2) 



(C2) 



wh ere the 7' s and their constituent functions are given 
in ( |3. 14 -3. 17). Here we shall numerically compute both 
factors as a function of impurity density and scattering 
potential. To facilitate the computation, it is convenient 



2Foln-^ ) l/i < 1 

^0, 



2Foln— ) 1^1 » 1 
^0, 



Born limit 
Unitary limit 



Thus, for a particular range of z, we shall set Y^caie to a 
value that allows us to evenly sample the transition from 
one limit to the other. In the end, these manipulations 
yield a new set of four dimensionless input parameters, 

{z,^, i?2,i?3}- 

The electrical and thermal/spin vertex corrections ob- 
tained for a typical set of input parameters, via the pro- 
cedure described above, are plotted in Fig. |^. Here we 
have assumed that T^(k) falls off slowly with increasing 
k (i?2 = 0.9, i?3 = 0.8) and have plotted the vertex 
corrections versus impurity fraction (z = 0.01% 1%) 
for seven scattering strengths from Born id — 0.001) to 
unitary [d = 0.999). Note that the electrical correction. 
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FIG. 5. Numerically calculated (a) electrical and (b) ther- 
mal/spin vertex corrections plotted as a function of impu- 
rity fraction (z = 0.01% — » 1%) for scattering strengths 
parametrized via (from bottom to top) d = {0.001 (Born), 
0.1, 0.3, 0.5, 0.7, 0.9, 0.999 (unitary)}. In all cases we have set 
Vacate = 20 and assumed that V(k) falls off slowly with |k 
(i?2 = 0.9, R3 = 0.8). 



f3yc — 1, can be quite significant while the thermal/spin 
correction, /Jj^;* — 1, is much smaller and vanishes as 
z — > 0. The difference between the two cases is seen 
most clearly in Fig. ^ where, for d — 0.001,0.5,0.999, we 
have replotted the electrical and thermal/spin correction 
factors on the same scale. On the scale of Pvc — 1, it is 
difficult to distinguish — 1 from the x-axis. Thus, 
we see that for all scattering strengths, the thermal/spin 
vertex correction is negligible compared to the electrical 
vertex correction. 

Additional insight is gained through consideration of 
the intermediate stages of the calculation which reveal 
that 



^(0) (0) (0) < 



1 



i 



(C6) 



In the Born limit, l/ln(po/ro) ^ z, which can truly be 
neglected for the small impurity fractions of interest. In 
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FIG. 6. Direct comparison of electrical {f3vc — 1) and ther- 
mal/spin (/Sjc" ~ 1) vertex corrections for (a) d = 0.001, (b) 
d = 0.5, and (c) d = 0.999. Note that in all cases, the ther- 
mal/spin correction is negligible on the scale of of the electri- 
cal correction. 



the unitary limit, l/ln(po/ro) ~ l/ln(l/z), which van- 
ishes much more slowly but is still small compared to 
terms of order unity. Thus, since we can treat the 7's as 
small quantities, the (mathematical) difference between 
the electrical and therm al/s pin cases is due to the manner 
in which the 7's enter ( |Cl| ) and (p2|). The 7's enter (3yc 
within order unity combinations of the form 7ln(po/ro) 
(and for some parameter values even 7ln^(po/ro)) which 
cannot be neglected. In contrast, the 7's enter /Jj^;* on 
their own (that is, in direct competition with terms of 
order unity) and can therefore be neglected for small z. 
Hence we say that electrical vertex corrections contribute 
even to zeroth order in the impurity density while ther- 
mal/spin vertex corrections do not. 



APPENDIX D: FERMI LIQUID CORRECTIONS 

In addition to the effects of vertex corrections, trans- 
port coefficients may be further modified due to under- 
lying Fermi liquid interactions between Landau quasi- 
particles (referred to hereafter as "electrons" to avoid 
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any confusion with the BogoHubov quasiparticle exci- 
tations of the superconductor). A detailed theory of 
the superfluid-J'erpii Uquid has been developed in the 
literatureJj'EJ'L2L3 In particular, a widely applicable ph&n 
nomenological approach has been devised by Leggett.Ea 
Within this approach, we consider three layers of Fermi 
liquid effects. (1) Mass renormalization: By virtue of 
our acknowledgment that "electrons" are really Landau 
quasiparticles, all masses should be viewed as effective 
masses, m* . (2) Current renormalization: Due to Fermi 
liquid interactions, the presence of a current yields an ad- 
ditional drag current resulting in an overall current renor- 
malization. (3) Response function modification: Fermi 
liquid interactions induce an effective "molecular field" 
which couples to the current-producing perturbation and 
modifies the response function, Kq{T), via 



Ko{T) ^ K{T) = 



1 + A/^o(r) 



(Dl) 



where A is a constant which depends on the interaction. 
In the case at hand, (1) has been built into the param- 
eters of our model, (2) yields a non-trivial effect which 
shall be considered below, and (3) can be neglected in the 
T — !■ limit with which we are concerned. This last state- 
ment follows because at low T, few quasiparticles are gen- 
erated, response functions (such as the normal fluid den- 
sity or conductiv ity) are small, and the higher order cor- 
rection terms in (Dl) are negligible. Thus, the dominant 
corrections to transport coefficients due to Fermi-liquid 
interactions enter simply via renormalization of the var- 
ious current density operators. To determine the nature 
of this current renormalization, we proceed as follows.Eil 
In the absence of Fermi liquid interactions, a generic 
(bare) current density operator takes the form 



(Note that in this context, the tilde denotes a dressed 
quantity, not a Nambu matrix). It is important to realize 
that the dispersion of every electron (not just the excited 
ones) is dressed as long as any electrons are excited. Al- 
though only the excited electrons interact, the dispersion 
measures the energy required to excite an electron and 
must therefore account for the interactions an electron 
would experience if it were excited. Also note that al- 
though Afe should also be modified due to Fermi Liquid 
interactions, we expect this effect to be less significant 
and shall therefore assume for the sake of simplicity that 
the gap is unaffected. Once the dispersion is dressed, we 
can Taylor expand in the change in to obtain a dressed 
equilibrium distribution function 



no(efc, Afc) = nQ{ek 



k' a.' 



(D5) 



Thus, the dressed current density operator takes the form 



where 



5ni 



(D6) 



(D7) 



is the difference between the true electron distribution 
and the equilibrium distribution in the presence of Fermi 
liquid interactions. 

In a superconductor, there are two types of elec- 
trons: those that compose the condensate of ground state 
pairs and those that form the Bogoliubov quasiparticles. 
Hence the equilibrium distribution function has both a 
condensate term and a quasiparticle term: 



jo = ^ gka^k Suka 

ka 



(D2) 



where is a velocity, gka is a coupling parameter (i.e. 
charge, spin, or energy), and 



no 



(D8) 



The equilibrium distribution of electrons in the conden- 
sate is just given by the coherence factor vj. so its deriva- 
tive is 



(D3) 



is the difference between the true electron distribution in 
the presence of the current- inducing perturbation, Uka, 
and the (bare) equilibrium distribution, no(£fci A^). (The 
bare designations refer to our neglect of Fermi liquid in- 
teractions aside from the extent to which they are in- 
cluded in the velocity via mass renormalization.) 

Once Fermi liquid interactions are turned on, excited 
electrons interact via the Landau interaction energy, 
f^l^ . In the presence of this interaction, the electron 
dispersion must be modified (or dressed) to account for 
the additional energy cost of interacting with other ex- 
cited electrons. 



fk = Cfe 



k' a.' 



fkk"'Sn 



k'a' 



(D4) 



5n£ 



d 



ifi-ii 

2 V Ek 



In the presence of impurities, the equilibrium quasiparti- 
cle distribution is given by the convolution of the quasi- 
particle spectral function (a Lorentzian of width Fq about 
-Efe) and the Fermi function np. Multiplying this by 
the probability the quasiparticle is an electron minus the 
probability it is a hole, (w^ — u^), we obtain the equilib- 
rium distribution of electrons that form the quasiparti- 
cles. In the T ^ Fq limit, its derivative is given by 



dr 



d 
dek 



To/tt 



Ek 
To/tt 



ElEl+Tl 



-El 



arctan 



np{Lo) 

El 

Ek 



(DIO) 
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Making use of these expressions and evaluating the 
dressed current density operator, 



Vfc 



Suka 



drin dn9 



dei, de 



(Dll) 



for the appropriate couphng parameters and velocities, 
the Fermi liquid renormalizations of the electrical, ther- 
mal, and spin currents can be computed. 
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